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. ABSTRACT 

Context. The VISTA Variable Survey (VVV) is able to map the Galaxy at / < 0° with an unpaired depth (at least 3 mag deeper than 
2MASS), opening new possibilities for studying the inner structure of the Milky Way. 
. Aims. In this paper we concentrate on the exploitation of these data to better understand the spatial disposition and distribution of the 

O structures present in the inner Milky Way, particularly the Long Bar and its interaction with the inner disc. 
Methods. To attain this, we present the Ks vs J-Ks diagrams obtained with VVV of regions with (-20° < / < -8°, \b\ < 2°). From 
them we derive the distribution of red clump giants with heliocentric distance as a proxy for the overall stellar structure of the Galaxy. 
Along with these diagrams, we analysed the distribution of photometrically selected red supergiants, in order to detect events of recent 
I ■ stellar formation. 

Results. The observations show the presence of a clear overdensity of stars with associated recent stellar formation that we interpret 
as the traces of the Long Bar, and we derive an angle for it of 41° ± 5° with line Sun-Galactic centre, touching the disc near / = 27° 
and I = -12°. The colour-magnitude diagrams presented here also show a lack of disc stars in several lines of sight, a fact that we 
associate with the truncation of the disc by the potential of this bar for R GC < 5 kpc. 

Key words. Stars: Supergiants - Infrared: stars - Galaxy: structure - Galaxy: stellar content 
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[ 1 . Introduction angle is around 25 deg. when low latitudes are excluded from the 

. fit, and around 45° in the plane regions within -15° < I < 30°. 

^ The signs that the Milky Way is a barred galaxy have increased Rather man two structures, iMartinez-Valnuesta & Gerhardl 
Q\ following the many discoveries about the morphology of the (|20TTh suggested that there is a scenario of a single struc- 
O . central regions of the Galaxy, most coming from detailed near- ture with a twisted major axis. We find this proposal accept- 
infrared (NIR) star counts. A bar-like distribution in the stars ab i e f rom a pure l y morphological point of view. However, bars 
7^ ; of our Galaxy is observed fro m the asymmetries in the infrared form as a result of an instability in differentially rotating discs 
> . surface brig htness maps (e. g. Ipwek et al.lll995l) and in source dSellwoodll 19811) . whereas bulges are a primordial galactic corn- 
counts (e.g. |Stanek et al.|[l994|;|Hammersley et al.|[1994j), which pon ent. Therefore, the scenario bulge+bar is not the same thing 
both systematically show more stars at positive Galactic longi- as a sing i e structure, and one should observe other differences 
tudes within < 30° and close to the Galactic plane compared to apal1 from the morphology. We think that the particular features 
the negative side. f this new proposal of a "single twisted bulge/bar" scenario 

leaves some observational facts unexplained, whereas the model 

While some papers refer to a thick triaxial structure (bulge f a misaligned bulge+long b ar successfully explains them all 

or thick bar) with a semi-major axis of around 2.5 kpc and dLonez-Corredoira et al.ll201lh . Furthermore, there are two dis- 

a position angle (PA) of 15°-30° with respect to the Sun- tinct populations in the thick bulge. The metal-rich population 

Galactic centre direction (e.g. [Ldpez-Corredoira et al. 2005; pre sents bar-like kinematics, while the metal-poor population 

Habi ngetalJ 1 20061: | Rattenbury et al. | | 2007fc | Vanhollebeke et al. | shows kinem atics corresponding to an old spheroid or a thick 

2009), other contributions suggest that there is also a disc, so the two main scenarios for the bulge formation co-exist 



X 



long thin bar, the in-plane bar, with 4 kpc in semi- wit hin the Milky Way bulge (Babusiau£itlD|Mi. 
major axis length and a posit i on angle of around 45° 
(e.g. | L6pez-Corredoira et aD 12007b ICabrera-Lavers et ai1l2007L We pay attention here to the long bar. We explore it in the 

Vall enari et al.ll2008t IChurchwell et al.ll2009l) . This long bar has negative longitudes (-20° < I < -8°, \b\ < 2°) with the data 

a tip in the positive longitude at the beg inning of the S cutum arm of 2MASS+VISTA-VVV in the NIR. For the positive longi- 

dDame et al.l fl986). As realized by ISevenster et al.1 (1 19991) . the tudes there are available studies including deep NIR surveys (e.g. 
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Cabrera- Lavers et al .I2007ll2008h . Four aspects are interesting to 
explore in this region: 

1. The stellar population of the bar itself, which can be 
followed by a tracer, such as the red clump population, 
as already carried out in the positive Galactic longitude s 



:ic longitudes 
1120071 20081. 



dHammerslev e t al. 20001 ICabrera-La vers et al. 

2. The tip of the long bar in the neg ative longitudes, which 
is ten tatively placed at / « -14° dLopez-Corredoira et alJ 
I2001h . Finding there a massive star formation region would 
be new supporting evidence for the presence of a long bar 
in our Gala xy, as it was for the tip of the bar in the positive 
longi t udes dLopez-Corredoira et al. 19991: Negueruela et alJ 
1201 U iGreen et all 1201 U iNguven Luong et al 1 1 1201 lb . Such 
regions can form because of the concentrations of shocked 
gas where a galactic bar meets a spiral arm, as observed at 
the ends of the bars of face-on external galaxies. 

3. The dust lane before the lon g bar in the negative galactic 
longitudes dCalbet et alJll996HMarshall et al.ll2008l) . 

4. The stellar density of the inner disc around and behind 
the long bar+bulge, and particularly the possibility that the 
in-plane regions have a deficit of stars, due to a trunca- 
tion in this density or to the effect of a flare in the disc 
dLopez-Corredoira et al.ll2004T) . 

All these features will be analysed using NIR colour-magnitude 
diagrams of the stellar populations, deriving the extinction, and 
separating within them the red clump stars (old population) and 
the red supergiants or very bright red giants (young population). 

This paper is linked to Amores et al. (2012), where the stellar 
counts produced with the same datasets used here in the form of 
CMDs are analysed. After a brief introduction of the datasets 
that will be used, a detailed description of the analysis tools we 
used follows in Sect. [3] In Sect. [4] the main results of the paper 
are given, leading to a summary of the more relevant conclusions 
in Sect. |5] 

2. Data 

2.1. VISTA-VVV 

The ESO public near-IR variability survey VVV is scanning the 
Milky Way bulge and an a djacent section of th e mid-plane with 
high star formation activity (Minn itTet al.l20l6T) . The survey will 
take 1929 hours of observations with the four metre VISTA tele- 
scope, across an area of approximately 550 square degrees in 
the inner galaxy. The final products will be a deep IR atlas in 
five passbands (ZYJHKs) and multi-epoch photometry in the 
TTs-band to build a catalogue of about a billion point sources, 
including about a million variable sources. The main goal of the 
VVV survey is to enable the construction of a 3-D map of the in- 
ner Milky Way using well-understood distance indicators, such 
as RR Lyrae, Cepheids and red clump giants. 

The first dataset, including observations from year 2010, has 
just been released. Sait o et all (120 1 2h discuss the observations, 
reductions, photometry, calibrations, and catalogues. For most 
regions of interest the limiting magnitude is m^ s ~ 18.0, and the 
image scale is 0.34 "/pixel. A more detailed description of the 
VVV data used can be found in Amores et al. (2012). 

2.2. 2MASS 

Although we are interested in the inner Milky Way, best viewed 
with VVV, we also study brighter populations, some of which 
might be too bright for this survey. We supplement the VVV data 



using 2MASS All-Sky Release dSkrutskie et al.l20 06). accessing 
its point source catalogue through the Gator interfac^J 



3. Data analysis 

3. 1 . Combining 2MASS and VVV 

We have made intensive use of the colour-magnitude diagram 
(CMD) instead of a source-by-source analysis, so a detailed 
combination of both datasets is not needed. Instead, for a more 
subtle approach we opted to adopt hard thresholds in K$ and /: 
above these thresholds we rely on 2MASS photometry, and be- 
low them VVV is the data source. 

Determining these thresholds requires careful calculation of 
the completeness of 2MASS. We restricted ourselves to -20° < 
/ < -8°, and in these regions completeness is dominated by 
source confusion, owing to the high stellar density. For each 
of the fields we built a histogram in magnitude, fit the brighter 
part with a second degree polynomial, and checked when the 
observed counts fall below 80% of the predicted value offered 
by this fit, adopting this value as our completeness limit. As we 
want to have a homogeneous set of CMDs, we selected a single 
set of thresholds for all our fields: (J, K s ) = (14, 12.5). 

As we are comparing data from two different surveys, with 
different instrumental configurations and calibration procedures, 
we also checked that both photometries are consistent. As both 
catalogues overlap in magnitude, for each field we crossmatched 
the common stars and checked for differences in magnitude. The 
typical values are around 0.05 mag in J and 0.02 mag in ^s, al- 
though in some fields subject to severe reddening these values 
can be higher, reaching differences up to 0.4 mag. In all our sub- 
sequent analysis, these systematic differences have been added 
to VVV magnitudes. 

We also checked the behaviour of the magnitude difference 
against colour and found no evident dependence, but in the red- 
der part of the CMD (with J - K$ > 4.5), there are few common 
stars, so this effect could be masked. 



3.2. The infrared colour magnitude diagram 

Since the advent of large-scale infrared surveys, the infrared 
CMD, usually of (J - K&) versus has become one of the 
standard tools for studying the structure of the Milky Way. 
Beyond the penetrating powers of this wavelength range that al- 
low us to reach the inner parts of our Galaxy, in these diagrams 
it is relatively easy to separate the two main visible popula- 
tions, dwarfs and giants. Although the morphology of the CMD 
changes greatly in different areas of the Galaxy, since they are 
heavily affected by interstellar extinction, in-plane fields offer a 
similar aspect, exemplified in Fig.Q] 

This CMD was obtained for a field centred at (/, b) = 
(-11°,0°). The main features that dominate the figure are out- 
lined in the right-hand panel: in blue and occupying the lower 
left-hand region of the plot, the dwarfs of the disc are the more 
abundant population, separated from the giants due to the extinc- 
tion and their intrinsic brightness. Even if caused by distance, 
dwarfs and giants could appear to us with the same apparent 
magnitudes, since giants are intrinsically brighter they would 
have to be farther away and be the subject of heavier interstellar 
extinction, since we are dealing with in-plane fields. This implies 
that they will appear shifted towards the red, occupying then the 
top right-hand part of the diagram. 



http://irsa.ipac.caltech.edu/applications/Gator/ 
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The effect of distance (that moves stars towards dimmer 
apparent magnitudes) and extinction (that makes stars appear 
dimmer and redder) results in an ideal population of stars with 
the same intrinsic magnitude and colour and evenly distributed 
along the line-of-sight that would appear as a diagonal line (not 
necessarily straight) in this CMD. The natural spread in mag- 
nitude of any real population changes this line into a strip of a 
certain width, and the uneven stellar density of the Galaxy trans- 
lates into overdensities in the (J - Ks,Ks) space. Because the 
later the type of a giant the brighter and redder it becomes, the 
lower envelope of this region is composed by late type G gi- 
ants, but the bulk of it will correspond to giants of spectral types 
around K2, also known as red clump giants (RCG), which are the 
most abundant giant population. The redder tail of Fig. Q] would 
be composed almost exclusively of M giants. The relative abun- 
dance of RCGs is high enough so that they appear as a distinct 
strip on the CMD while later type giants, much less abundant, 
are visible towards redder and brighter magnitudes. 

To build the CMDs used in this work (including the one from 
Fig. [I), we always use 0.5° x 0.5° fields centered on the nominal 
coordinates. As the source density is high, instead of plotting 
star by star, we built a density diagram using tiles of 0.1 mag 
in colour and magnitude, as can be seen in Fig. Q] Also, the 
variation range of N, the number of stars in each one of these 
tiles, is wide, so that to appreciate all the substructure present 
in these CMD it is mandatory to plot instead some function that 
compresses this range. After a lengthy trial and error process we 
opted for the function 



/ = 300 x log (N + l) + N, 



(1) 



where N is the number of stars per tile and / the plotted intensity. 
This selection combines the compression factor of a logarithmic 
intensity plot but without losing the ability to show density vari- 
ations of ~ 100 stars per tile, changes that in a pure logarithmic 
plot would be invisible. As our analysis of these density vari- 
ations was purely qualitative, we opted to use a function that 
maximizes the visible information instead of one with immedi- 
ate physical translation, such as a square root or a logarithm. 

The resulting CMDs, used in subsequent sections, can be 
seen in AppendixlAl 

3.3. The red clump method 

As said above, the ubiquity of the RCGs makes the so-called 
red clump method one of the standard tools for analysing 
CMDs i n the infrared. Its use a n d ap plication is discussed in 
depth in Lope z-Corredoira et al.l d2002l) (hereafter LC02), but 
in this section we offer a brief outline of its application to the 
2MASS+VVV composite CMDs, as the results of this method 
is used thoroughly through the text. 

These RCGs have a very narrow luminosity function that it 
is also calibrated well (LC02), so their intrinsic colour and lu- 
minosity are well known. It is possible then to transform their 
position in the (/ - Ks,Ks) plane to (iJ,Ak), once we assume 
an extinction law (in the form of the relative extinction coef- 
ficient Aj/Ak). A more detailed d escription of the method and 
the errors involved can be found in ICabrera-Lavers et alJ (T2008) 
and references therein. Following this study, we use (J - Ks)o = 
0.70 ± 0.05, M K = -1.62 + 0.03 andAj/A K = 2.518. 

To perform this calculation there is a previous step that needs 
to be taken: the calculation of the (J - Ks, Ks) curve that repre- 
sents the distribution of RCGs over the CMD. The simplest way 
to do this is to cut the CMD into several strips in Ks and then fit 



the (J-Ks) distribution of each strip with the sum of a second de- 
gree polynomial and a Gaussian; the /i of this gaussian gives the 
fiducial point that represents the locus of the RCG for that strip 
(LC02). The problem here is that, in the presence of high inter- 
stellar extinction (see for example the CMDs for b — 0° in Fig. 
IA.4t , the RCG strip almost appears horizontal, and so projecting 
the data of each strip over the (J - Ks) is not optimal, because it 
mixes stars with very different colours but with roughly the same 
Ks magnitude, resulting in a loss of information. 

To overcome this, as a first step in our analysis, a rough fit of 
the RCG strip is performed by eye, and this (J - Ks,Ks)o curve 
is used as the entry point of an automated process that, instead 
of performing cuts over Ks and projecting them over the {J-Ks) 
axis, defines a local coordinate system (X, Y) that is perpendicu- 
lar to this (J - Ks,Ks)o curve at each point. As these new axes 
take the reddening vector into account, stars with similar magni- 
tudes but different reddening appear separated when projecting 
over either X or Y, so a much more refined (J-Ks, Ks) curve 
can be derived. 

Once we have the variation in (J-Ks) with Ks for the RCGs, 
we can extract this population from the CMD. Although this can 
be optimized for each line-of-sight, and because our study re- 
quires the comparison of data at different Galactic longitudes, 
we consider that everything within 0.5 mag in colour and mag- 
nitude from this curve is a RCG. 

3.4. Photometric search for massive stars 

One of the best tracers of stellar formation regions are red su- 
pergiants (RSG). Since these stars are very bright in the in frared 
(with a conservative limit of Mk < -8, IWood et al.l lT983). they 
shine through the obscuration present in the plane. One of the 
expected effects of a galactic bar is a triggering of star formation 
at its tips and, depe nding on the bar strength, over its dust lanes 
(see, for example. ISheth et ani2002l) . As such, the study of the 
distribution of RSGs with Galactic longitude could give us clues 
on the structure of the bar. 

In the near end of the bar, around I ~ 28 ° several clusters 
of these stars have already been detected (i.e. [N egueruela et all 
1201 It [Gonzalez-Fernandez & Negueruela 2012). In these clus- 
ters, RSGs will have apparent magnitudes 7 < < 4. If we as- 
sume that the bar forms a 45° angle with the Sun-Galactic centr e 
line and has a 4 kpc semimajor axis (ICabrera-Lavers et al.l2 008). 
this would put the far end of the bar at ~ 1 1 .2 kpc, while the near 
one sits at ~ 5.9 kpc. This difference, in distance modulus, is of 
~ 1 .4 mag, so the same RSG that we see at the near end of the 
bar would have magnitudes in the range 9 < < 6, allowing 
for a ~ 0.5 mag in extra extinction due to the presence of the 
dust lane of the bar. 

It is worth noting that, although we have selected a <p — 45° 
bar, ot her possible position angles have been postulated (for ex- 
ample, |Babusiaux^Gito^|2005l with $ = 22° + 5.5°), but 
always with <p < 45°. This implies that for these bar configu- 
rations, the far end should be farther away than for a 45° angle 
bar. The difference in brightness for a given RSG between both 
configurations is only of 0.1 mag (corresponding to 0.6 kpc), 
and our selection criteria for RSG is therefore not biased against 
these models. 

For this reason, the possible RSGs situated at the far end 
of the bar would be very difficult to separate from the popu- 
lation of late type giants (they also have very similar intrinsic 
colours) of the disc. One possible way to separate this popu- 
lation is to use the reddening fre e pseudo-colour Q, defined a s 
Q = (J - H) - 1.8 x (H - Ks) (Negueruela & Schurd3 l2007h . 
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Fig.l. Left: Sample CMD of a 0.5° x 0.5° field around (l,b) = (-11°, 0°), combining data from VVV and 2MASS. The white 
dashed line marks the transition between both surveys. Right: Sketch of the same plot with the regions that occupy dwarfs (blue), 
RCGs (yellow), and later-type giants (red), as discussed in the text. 



Generally speaking, the later (and intrinsically redder) a star is, 
the higher value of Q it will have, with late type giants (types 
beyond K0) having Q > 0.41. RSGs, because their atmospheres 
are very extended, complex, and often present mass loss, will 
not deredden properly, and so they roughly occupy the range 
0.1 < Q < 0.4 dMessineo et al.ll201 

Another applicable threshold affects the colour of the candi- 
dates. RSGs are intrinsically red sources, with intrinsic colours 
(J - Ks)o > 0.7. Since we are also looking at sources under 
heavy extinction (as is expected for the inner Milky Way), we 
can safely assume that any RSG will have (J - K&) > 1.3. This 
step eliminates late type dwarfs from the disc that will fulfill 
most of the previous criteria but, as they are located on the local 
disc will have lower colour excesses. 

Summing up, the criteria that we used to select RSG candi- 
dates over 2MASS data are 

1. 0.1 < Q < 0.4, 

2. m K < 9, 

3. (J-K S )>1.3. 

Negueruel a et al.l d201 ll) performed a spectroscopic survey 
around the massive cluster RSGC3, confirming the nature of 
stars selected based on these criteria, and find that the initial sam- 
ple is entirely composed of RSGs and very evolved giants, with 
stellar types later than M5. A 38% fraction of the objects turned 
out to be class I sources (supergiants), while 28% were class II 
(bright giants), and the remaining 34% identified as late class III 
giants. 




-8 -10 -12 -14 -16 -18 
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Fig. 3. Number of red giants per square degree as a function of 
Galactic longitude and distance for fields with b = 0°. To in- 
crease contrast, the plotted quantity is the square root of this 
number, N 1 ^ 2 . The white dashed line represents the complete- 
ness limit at each line-of-sight. The blue and red dashed lines 
mark the place of a 45° bar and a ring of r — 4.5 kpc, respec- 
tively. The sudden drop at I - -8° is due to the combination 
of higher reddening and stellar density. All these giants verify 
|fc|<0.5°. 



4. Results 

4.1. Structure of the Galactic bar 

One of the most striking facts when analysing the CMDs pro- 
duced with 2MASS+VVV for fields with |Z| < 0° and b = 0° 
is that there is a sudden drop in the number of giants around 
I = -13° (see Fig. |2|. Since the total number of sources per 
square degree would be dominated by the dwarf population of 
the disc and we are more interested in the innermost region of 
the Milky Way, we focus our study on the giants, that reach be- 
yond the disc within the limiting magnitude of VVV. 

Using the method described in Sect. 13.31 we can extract the 
RCGs from a CMD and estimate their distance. This allows us 



to reconstruct the in-plane 2-D distribution of these stars, as can 
be seen in fig. [3] Beyond I = -8°, where extinction and source 
density affect our sample, the behaviour of the number of giants 
with (I, d) is rather smooth, showing traces of a rather straight 
structure that ends beyond I ~ -12°, where there is an appre- 
ciable drop in the number of visible giants. This could easily be 
associated with a bar-like structure, but before doing this there 
are some issues that need to be tackled. This drop in the number 
of giants visible in Fig.|2]could be due to an increment in the ex- 
tinction along the line-of-sight, but this is unlikely: in this same 
figure it can be seen that the morphology of the RCG strip is sim- 
ilar in all the lines of sight, with higher latitudes and lower ex- 
tinction. This implies a similar (d, Ak) behaviour for all of them. 
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Fig. 2. CMDs of the in plane fields (b = 0°) between I - -12.0° and I = -13.5°. The colour scale is the same as Fig.Q~| and the white 
dashed line represents the transition from 2MASS to VVV data. The rather sudden drop in the number of stars between I = -12.5° 
and I = -13.0° is evident in this plot. 



Using the method described in Sect. 13.31 we can estimate the 
variation of E(J - Ks) with d, plotted in Fig. [4] Comparing it 
with Fig. |2] we see that, although the extinction along the line- 
of-sight is almost equal for both of them, the number of visible 
giants is appreciably lower for I = -13°. 

The density distribution of Fig.[3]could also be an effect of 
the variations in the limiting magnitude with Galactic longitude. 
A simple exponential density distribution, when observed with 
the instrumental setup of VVV, would show a maximum in the 
number of stars near completeness, for when we approach this 
limit we start to gradually lose more visible stars, modulating 
the subjacent stellar distribution. As hinted before, in off-plane 
fields these effects are less severe, and in Fig.|5]we can see how 
this bar-like structure is visible, well before the completeness 
limit is reached. 

This structure corresponds very well to the exten- 
sion of the Long B ar, as has been observ ed previously 
in inte grated light (iHammerslev et all l2000h and stellar 
counts (Lopez-Corredoi ra et al.l 120011) . to the fourth quadrant. 
Assuming that this bar is symmetric with respect to the Galactic 
centre, the positions of its tips provide a strong constraint on 
its size and orientation. We place its far end at / = -13° (as 
suggested by Figs. [3] |9l and [5), and the nearby end in / = 27° 
dCabrera-Lavers et al.ll2008l) . Although these measurements are 




0.0 Lj , , , i , , , i , , , i , , , J 

2 4 6 8 10 

d (kpc) 

Fig. 4. Variation in E(J - K s ) with the distance along the line- 
of-sight for the in-plane fields at I = -12° (black) and I = -13° 
(red). The strange behaviour seen in this last curve at ~ 4.5 kpc 
is due to the imperfect combination of 2MASS and VVV data 
in this line-of-sight. For an explanation of the negative slope in 
parts of the curves, see Sect. 14.31 
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not very precise, the error should not exceed +1°. Assuming 
that R C c = 8.0 ± 0.5 kpc, we obtain an angle of 41° ±5° and 
a total len gth of 8.2 + 0.6 kpc. This is c l ose to the values cal- 
culated by Lopez-Corredoira et alJ d200ll) : ICabrera-Lavers et al.l 
(120071 2008). but clearly differe nt from the proposed geometry 
of Babusiaux & Gilmore (2005). 

The structure in Figs. [3] and [5] seems to have a wider angle, 
close to 45° but the range in galactic longitudes is small, and 
these differences in orientation are difficult to judge. Also, the 
derivation of d for this figure is affected by extinction, which 
could warp the visible structure if there is variation between dif- 
ferent lines of sight. The properties derived using only the angu- 
lar positions of the bar ends are very robust, not affected by any 
of these effects, the distance to the Galactic centre being the only 
assumption made. 

Another interesting feature that, for example, can be ob- 
served in the CMDs of Fig. [2] and the extinction curves from 
Fig. [4] is that there is a huge amount of interstellar material con- 
centrated in a narrow strip around 9 kpc at I = -13° (although 
this distance estimate varies with the assumed extinction law, see 
Sect. 14.3b . This wall of extinction makes the RCGs appear as an 
almost horizontal strip on the CMD, but its effect is constrained 
to the plane. In Fig. [6] we plot the diagrams of a vertical cut at 
/ = -12.5° (near the end of the bar), and we can see clearly how 
the effect of interstellar extinction decreases from b = 0.0° to 
b = -0.5°. This feature of the interstellar medium can be clearly 
followed in the plane from I — -13° to / = -8° (see Figs. lA.ll to 
lA~T2l . 

In these longitudes it is easy to see that this structure is too 
far away to be mistaken for the 3-kpc arm/molecular ring (con- 
trary to what might happen in the first quadrant), so this sudden 
increase in the extinction (hence, the density of interstellar ma- 
terial) is a clear clue to the existence of a leading dust lane in 
the Long Bar of the Milky Way. A similar structure was pre- 
viously detected in the inner Milky Way ov er extinction maps 
and molecular emission dMarshall et al.l 2008). although at lower 
|/|. Dust lanes are a co mmon companion to stellar bars in other 
galaxies (for example, IComeron et al. 1 120091) . so the presence of 
such a structure in our Galaxy is a hint of the existence of a bar 
up until / ~ -12°. 

As we move outwards, this dust lane disappears, but the 
interstellar extinction grows more or less continuosly for I < 
-14.5°. As can be seen in Fig. [7] the variation in E(J - K$) 
with d is quite different from that of the dust lane. Instead of a 
sudden increase in the amount of interstellar material, the vari- 
ation of Fig. [7] is much smoother. This matches very well what 
is expected of a near-tangent line-of-sight to a gas rich structure 
such as an arm (see Sect. 14 . 4b - 



4.2. Structure of the inner disc 

A careful look at Fig.[5]presents us with another interesting fea- 
ture: the depletion of giants present from -18° < / < -12° 
(particularly evident at b = 1.5°). This truncation of the stel- 
lar density of the disc is visible in every galactic latitude ex- 
cept for \b\ < 0.5°, due to completeness (see Fig. [3]). This is 
one of the strongest visible eff ects of a bar; its typical poten- 
tial (Binnev & Tremainej 120081) has two maxima in the direc- 
tion of its min or axis (labelled L4 and L 5 in Fig. [8] following 
the notation of lBinnev & fr emaine 2008). These maxima sweep 
away stars and gas outwards, effectively emptying the disc in- 
side the corotation radius. In Fig. [5] this appears clearly in the 
first plot, and there are traces of the same effect in the third, 
but masked mostly by the effects of extinction. The method de- 



scribed in Sect. [373] assumes that the extinction varies smoothly 
along the line-of-sight, so the dumpiness of the real distribution 
of Ak with d translates into a cumulative blur of the features in 
Figs.[3|and|5] 

This truncation is not only visible in the stellar content, but 
it also affects interstellar matter. As can be seen in FigJU from 
~ 4.5 to ~ 9 kpc there is little variation in E(J - K$) with d. As 
the colour excess is proportional to Ak s and this quantity, in turn, 
is a function of the integral of the volume of interstellar matter 
along the line-of-sight, a flat curve implies that there is little to 
no underlying gas, matching the behaviour of the giants in Fig|5] 

Even so, the photometric depth of VVV is high enough to 
probe the far disc for the first time. In some of the CMDs of 
appendix lAl the trace of the RCG runs vertically for tok > 12 
(particularly evident for (I, b) = (-12.5, -1 .0)). This trace shows 
a secondary maximum, visible on several lines of sight (as, for 
example, in Fig.|6]l. These maxima appear around 13 kpc, and are 
clearly visible in the reconstructed distribution of Fig. [5] (more 
evident in the middle and bottom panels). 

This increase in stellar density shows us where the disc re- 
gains its density beyond coro tation. This effect has also been 
noted in lGonzalez et alJ (1201 ll) at \b\ - 1°, although the authors 
interpret these maxima as traces for an arm at 1 1 .2 kpc. This 
seems unlikely, since a one-degree separation from the plane 
implies a height of 200 pc at 1 1 kpc (300 pc for a 1.5° line-of- 
sight), and arms are structures enclosed in the thin disc, so have 
very small height-scales. On top of this, the effective volume of 
an arm that a line-of-sight crossing it probes is very small unless 
at a tangent point. Since these structures are mostly composed of 
young populations, they are not traced very well with RCGs (a 
late population, with r ~ 2 Gyr), and it is very unlikely that they 
could appear clearly at any line-of-sight, as in Fig. [5] 

4.3. The behaviour of interstellar extinction 

Many of the previous results depend, or at least are affected by, 
extinction. In the derivation of Figs. [3] and [5] we have assumed 
a simple extinction law that remains constant with d, I, and b 
and that yields Aj/Ak = 2.518. As can be seen in Fig. [4] beyond 
9 kpc, the slope of the (d, E(J - Ks)) curve becomes negative. 
This is physically impossible, and is a direct effect of this choice 
of extinction law, for the derivation of d depends on the calcula- 
tion of Ak and, as follows: 

E(J - tf s ) = (J - K S ) -(J- K s ) 



E( J - K S ) — Ai — Ak s = Ak s -7— — 1 

This implies: 

= (J - K S ) -(J- K s ) 
Ks Aj/A Ks - 1 

Once we isolate a given population of stars that we will use 
as standard candles, from two measurables, (J-Ks) and wk s , we 
therefore derive d, E(J - Ks), and Ak s - To complete the recipe 
we need to know, apart from the aforementioned Aj/Ak s , the 
luminosity function of the population under study, through its 
{J - Ks)o and Mk s . While a variation in any of these three con- 
stants would change Fig. [4] we know that: 

- An offset in (J - Ks)o or Mk s would, in exchange, offset the 
curves in Fig. |4j but not change their slope or overall shape. 

- For the case of the RCGs, there is a known dependen cy of 
(J - K s )o with metallicity (ICabrera-Lavers et al]l2007l) . but 
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we can see in Fig. [4] that, since the extinction increases in 
such a quick manner, most of the variation (particularly the 
part in which the curve acquires negative slope) takes place 
in less than a kiloparsec, so the metallicity spread is expected 
to be low. For this reason, while the overall shape of the d vs. 
E{ J-Ks) curve might change due to this effect, this variation 
would affect the slope of the nearby part of the graph with 
respect to the points that are more distant from the observer, 
but it will not change the slope measured for d > 8 kpc. 

Taking these points into consideration, the only remaining 
possibility to explain this negative slope is our choice of a value 
of Aj /Ak s not representative of the observed physical conditions. 
We can tinker with this value in order to bring the curves in Fig. 
[4] into coherence with the expected behaviour of extinction with 
distance (that, by definition, needs to be a monotonically increas- 
ing function), obtaining a lower limit for the observed value of 
this ratio. We find that 

A> 2 .9, 

so far away from the ty pical values of iRieke & Lebofskvl d 19851) 
or lCardelli et all (1 19891) . both close to Aj/A Ks = 2.5. This points 
to an interstellar material that is more transparent towards the 
red than predicted by the extincti on law, more in concord ance 
with the results from studies as Nishiyama et al. (2009) and 
IStead & Hoarel (120091) . This effect is not exclusive to the lines 
of sight used to produce Fig. |4j heavy extinction transforms the 
RCG strip in a CMD into a nearly horizontal band, so similar 
behaviors are observed in all the fields with -8.0° < I < -13° 
and b = 0.0° in Figs.lXHto [Q2l 

A change in the Aj/Ak s extinction ratio implies an expan- 
sion/contraction in the distance scale in Figs. [3] and [5] but will 
not change the overall shape of the structural features discussed. 
Because the analysis of these figures and the results drawn from 
are qualitative and the physical characteristics of the bar are cal- 
culated from the Galactic longitudes of its ends, a measure that 
is independent of extinction, our results are not affected by the 
choice of extinction law. Since a proper determination of the 
value of Aj /Ak s is beyond the scope of this w ork, we opt to re- 
tain the value from Rieke & Lebofskvl (1 19851) for homogene ity 
with similar studies (for example. ICabrera-Lavers et al.ll2007l) 

This distance scale factor introduced through Aj /Ak s can be 
appreciated when comparing Fig. [3] with Fig. [5] the bar over- 
density seems to shift away from the observer at higher \b\. As 
increasing Aj/Ak s implies decreasing Ak s and a higher derived 
d, the error in this last parameter due to the difference in Aj/Ak s 
will be a monotonically increasing function of E(J - Ks), hence 
higher for in-plane fields. 

4.4. Bar-driven stellar formation in the Milky Way 

A structure such as the Long Bar is expected to affect its envi- 
ronment in several ways. One of them is the triggering of star 
formation where its tips interact with the Galactic disc. Since 
RSGs are massive stars, their lifespans are expected to be short, 
just a few tens of Myr, so they trace recent stellar formation. We 
built a sample of photometrically selected RSG candidates fol- 
lowing the methods discussed in Sect. 13.41 namely, we selected 
stars with 0.1 < Q < 0.4, ( J - K s ) > 1.3 and m K < 9. Since 
these values are well within the completeness of 2MASS, we 
use this survey. Doing so lets us obtain candidates for all the 
inner Galaxy, from I = 40° to / = -30°. 



The vertical distribution of RSG candidates is, as expected 
for a young population, concentrated in low b fields, and it can 
be roughly approximated with a Gaussian of <x = 2.4°. We find 
that in the region with -30° < / < 40° there are close to 30000 
candidates with \b\ < 1° (roughly a density of 200 deg~ 2 ) and 
60 000 in off-plane fields. 

As reasoned in ICabrera-Lavers et al. (2007), the long bar is 
a very flat structure, so to locate traces of its associated stellar 
formation it is natural to restrict ourselves to the Galactic plane. 
Since several structures coexist in the region with \l\ < 30°, we 
compared the in-plane (\b\ < 1°) distribution of candidate RSGs 
with off-plane (1° < \b\ < 90°) fields. The results can be seen 
in Fig.|9](to facilitate comparison, the off-plane distribution has 
been normalized to have the same total number of objects as for 
in-plane fields). As can be seen, rather than being smooth, the in- 
plane distribution presents a series of overdensities not present in 
off-plane fields (or present in a much less conspicuous manner). 

These overdensities (outside the Galactic center) give us 
hints to possible underlying structures. This distribution follows 
the 2.2 /mi and 12 fim emission studied by Hammerslev et al. 
(1 19941) closely. As the structure of the first quadrant of the 
Galaxy is much better studied, the nature of some of these peaks 
has already been established: the overdensity at I = 32° is re- 
lated to the tangent line-of-sight towards the Scutum-Crux arm, 
and its c ounterpart in the fourth quadrant was already estab- 
lished by iHammerslev et al.l (1 19941) beyond I < -29° (hence 
outside our analysis). The next peak at I = 28° is known to 
be related to one of the more massive episodes of star forma- 
tion in our Galaxy, containing several massive clusters with 
M > 10 4 Mq, an d probably with a combined mass well 
in excess of 10 s M n (|Gonzalez-Fernandez & Neguerue ial l2(Hl 
Negueruela et al. 20 fH). Since the wings of this peak are wide, it 
also contains the extremely ma ssive clusters RSGC 1 (|Figer et aL 
2006; iNegueruela et al.l 12010) and Stephenson 2 dDavies et al. 
2007). All these clusters have been related to the interaction of 
the Long Bar and the Scutum-Crux Arm/molecular ring; this 
area also contains the more massive star form ing region of the 
Milky Way, W43 dNguven Luong et alJ[201 ll) (although it does 
not harbour any RSG yet). 

It is expected that a similar wealth of structures should occur 
at the far end of the bar, here identified at I = -13°, although 
the associated stellar formation has yet to be discovered. Some 
recent studies of massive clusters have started this search, but so 
far the only confirm ed massive association detected is Mercer 8 1 
dDavies et alJl2012h . It is associated by the authors to the far tip 
of the bar, but sitting at I = -21.6° is likely to be part of the 3- 
kpc arm, with a tangent at / = -22° (LC01). This is a prominent 
feature in gas e mission that was the first indication of our Galaxy 
being barred (Ide Vaucouleursl 19641) . Although it was previously 
suspecte d that this arm has not bee n subject to recent stellar for- 
mation (Dame & Thaddeus 2008), in Fig. [9] it appears clearly 
delineated through its tangent line-of-sight, pointing towards an 
associated RSG population. 

This arm occupies roughly the same space as the so-called 
molecular ring, a fe ature visible mostly in molecular emission 
dJackson et al.l 2006). Both structures are mostly young and will 
be poorly traced by RCGs. As a result, only their tangents are 
visible in our data and, it is visible in Fig. [5] is very difficult 
to distinguish between the two possible configurations, a closed 
ring or two tightly wound arms, by looking at both our CMDs 
and RSG distribution. 

The exact stellar nature of the remaining overdensity in the 
near end, visible in Fig.|9]spanning from I ~ 20° to / ~ 15°, has 
not been established yet. It might be part of the Scutum-Crux 
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arm, might be stellar formation of the bar itself B or a part of the 
molecular ring or the 3-kpc arm( Hammersle y" st al.ll994h (or any 
combination of the three). It does not have a clear counterpart at 
the far end, and this could point to it being part of the 3-kpc arm. 
As can be seen in Fig. [8] the cross section offered by the far pass 
of the arm to the line-of-sight is much wider, so all the stellar 
formation that in the first quadrant is visible in 20° < Z < 10° 
would appear, due to the effect of projection, around the tan- 
gential line-of-sight towards the arm in the fourth quadrant (at 
I = -22°). If there is any stellar formation along the dust lane, 
in Fig. [9] it would be visible with / < -13° and so it would be 
mixed with that of the tip of the Long Bar and the increased den- 
sity towards the Galactic centre. As such, it is almost impossible 
to distinguish between both hypotheses without a proper spec- 
troscopic study of the individual sources themselves. 

5. Conclusions 

In this work we obtained combined CMDs of 2MASS and VVV, 
mapping for the first time the innermost Milky Way at negative 
galactic longitudes. Using these data along with photometrically 
selected RSGs -used as proxies for recent stellar formation- we 
mapped the structures present at these latitudes, offering a fresh 
picture of the structure of our Galaxy. 
The main conclusions drawn follow: 

- Based on the analysis of the CMDs and the 2-D distribution 
of the RCGs along several lines of sight, we showed there 
is a flat structure that runs up to / ~ -12°. It corresponds to 
the extension of the Long Bar into the fourth quadrant. With 
this new constraint, we concluded that it forms an angle of 
41° + 5° with the Sun-Galactic centre line. 

- This Long Bar affects the structure of the disc. One of the 
more severe effects is the truncation of the stellar density 
with Rqc - 5 kpc. Because it corresponds with two max- 
ima of the bar potential, this area is devoid of stars and gas 
(outside the body of the bar itself). 

- Using photometrically selected RSGs as a proxy for recent 
stellar formation, we demonstrated another of the predicted 
effects of a galactic bar, the triggering of formation regions at 

its ends, clearly visible at I ~ 27° and I 12°. The tangents 

to the Scutum-Crux and the so-called 3-kpc arm are also vis- 
ible in the distribution of these stars, showing that in both 
arms there are traces of recent massive stellar formation. 



Appendix A: VVV+2MASS CMDs of the inner Milky 
Way 
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Ftig. A.2. CMDs of the in plane fields with b < 0°. The 



axis and colour scale 
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Fig. A.10. CMDs of the in plane fields with b > 0°. The axis and colour scale are those of Fig.Q] 
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Fig. 5. Same plot as Fig. [3j but for RCGs with -0.75° > b > 
-1.25°, (bottom), 0.75° < b < 1.25° (middle), and 1.25° < b < 
1.75° (top). Notice the difference in scales within these plots and 
also with Fig. [3] 
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Fig. 6. CMDs of the fields with / = -12.0° between b = -1.0° 
and b = 0.0°. The colour scale is the same as Fig.Q] As can be 
seen, there is a sudden increase in the extinction at (J-Ks,m^) ~ 
(2, 14) in the plane that is not visible at higher \b\. Also, the sec- 
ond maxima of the RCGs is visible in the bottom plot, roughly 
at (J-K 5 ,K S ) = (2, 15) 




Fig. 7. Same as Fig.|4]but for the in-plane fields at I = -18° and 
/ = -19°. 
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Fig. 8. Schematic view of the inner Milky Way, showing two 
possible configurations: a ring (red dashed line) and two arms 
(solid black lines). While both are almost indistinguishable in 
the first quadrant, it is possible to distinguish between them in 
the far end of the bar and the tangent line to the arm. The light- 
red shaded areas represent the areas where the maxima of the 
potential of the bar (L4 and L5) sweep away stars and gas. The 
bulge has been left out of this representation. 
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Fig. 9. Number of candidate RSGs detected in 2MASS, both for 
fields in the Galactic plane Qb\ < 1°) and outside (1° < \b\ < 
90°). To enable easy comparison, both measurements have been 
normalized to the same total number of stars. 
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